Compressible Turbulence: The Cascade and its Locality 



Hussein Aluie 

Applied Mathematics and Plasma Physics (T-5) & Center for Non-linear Studies, 
Los Alamos National Laboratory, MS-B258 Los Alamos, NM 87545, USA 
(Dated: January 4, 2011) 

We prove that inter-scale transfer of kinetic energy in compressible turbulence is dominated by 
local interactions. In particular, our results preclude direct transfer of kinetic energy from large- 
scales directly to dissipation scales, such as into shocks, in high Reynolds number turbulence as 
is commonly believed. Our assumptions on the scaling of structure functions are weak and enjoy 
compelling empirical support. Under a stronger assumption on pressure dilatation co-spectrum, 
we show that mean kinetic and internal energy budgets statistically decouple beyond a transitional 
"conversion" range. Our analysis establishes the existence of an ensuing inertial range over which 
mean SGS kinetic energy flux becomes constant, independent of scale. Over this inertial range, 
mean kinetic energy cascades locally and in a conservative fashion, despite not being an invariant. 

PACS numbers: 47.27.-i, 47.40.-x, 47.27.E-, 47.27.Jv, 47.27.em 



Turbulence is a phenomenon that pervades most liq- 
uid, gas, and plasma flows in engineering and nature, 
ranging from high-speed engines, nuclear fusion power 
reactors, and spacecraft re-entry, to star formation in 
molecular clouds, and supernovae. While the traditional 
Richardson-Kolmogorov-Onsager picture is a successful 
theory of incompressible turbulence, all aforementioned 
systems are characterized by significant compressibility 
effects. We have recently laid a rigorous framework in [1| 
to study scale-coupling in compressible flows and to an- 
alyze transfer of kinetic energy between different scales. 
The purpose of this Letter is to explore if such trans- 
fer of energy takes place through a cascade process and 
whether the cascade is scale-local. 

Kolmogorov's 1941 theory of incompressible turbu- 
lence makes the fundamental assumption of a scale-local 
cascade process in which modes all of a comparable scale 
^ £ participate predominantly in the transfer of energy 
across scale £. If, furthermore, the cascade steps are 
chaotic processes then it is expected that any "memory" 
of large-scale particulars of the system, such as geometry 
and large-scale statistics, or the specifics of microscopic 
dissipation, will be "forgotten." This gives rise to an in- 
ertial scale-range over which turbulent fluctuations have 
universal statistics and the flow evolves under its own in- 
ternal dynamics without direct communication with the 
largest or smallest scales in the system. 

Therefore, scale locality of the cascade is crucial to 
justify the existence of universal statistics and to war- 
rant the concept of an inertial range. It is, furthermore, 
necessary for the physical foundation of large-eddy sim- 
ulation (LES) modelling of turbulence. It motivates the 
belief that models of subscale terms in the equations for 
large-scales can be of general utility, independent of the 
particulars of turbulent flows under study. While scale- 
locality in incompressible turbulence stands on firm the- 
oretical [1, [3l and numerical 0, 0] grounds, no similar 
results exist for compressible turbulence. In fact, there is 



a widespread belief especially common in the astrophys- 
ical literature which maintains that a "finite portion" of 
energy at a given scale must be dissipated directly into 
shocks through non-local transfer in scale (see for exam- 
ple Q). Moreover, the idea of a cascade itself is without 
physical basis since kinetic energy is not a global invari- 
ant of the inviscid dynamics. Hence, the notion of an in- 
ertial cascade-range in compressible turbulence remains 
tenuous and unsubstantiated. 

In this Letter, we prove under modest assumptions 
that transfer of kinetic energy is indeed local in scale. 
Under a stronger assumption, we will further show that 
kinetic energy cascades conservatively despite not being 
an invariant. We reach these results by a direct analysis 
of the compressible Navier Stokes equations, without use 
of any closure approximation. The equations are those 
of continuity and momentum: 

dtp + V-{pu}^0, (1) 

dt{pu) + V-(/9uu) = - VP + /iV-(Vu + i V-uI) + pi, (2) 

o 

and either internal or total energy, supplemented with 
an equation of state for the fluid. Here, f is an external 
acceleration field stirring the fluid, and we have assumed 
a constant dynamic viscosity, p,. 

Our analysis is based on a coarse-graining (or filter- 
ing) approach expounded in [l|, in which we observed 
that any scale decomposition aimed at studying inertial- 
range dynamics must satisfy an inviscid criterion, i.e. 
it must guarantee that viscous momentum diffusion and 
kinetic energy dissipation are negligible at large-scales. 
We proved that a Favre decomposition meets such a 
requirement. Using classically filtered fields, a£(x) = 
Jd^r Gf(r)a(x-Hr), with kernel Gi{r) = e-^G{r/e) that 
is smooth and decays sufficiently rapidly for large r, a 
Favre filtered field is weighted by density as af(x) = 
pa£(x)/p£(x). From coarse-grained continuity and mo- 
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mentum equations, wc can write down a large-scale ki- 
netic energy budget, 

+ V-J, = -He - A, + P.V-u, -Di + e"^". (3) 

Eq. ([3]) describes instantaneous kinetic energy evolution 
at every point x in the flow and at scales > £, for ar- 
bitrary £. Our approach, therefore, allows the simulta- 
neous resolution of dynamics both in scale and in space, 
and admits intuitive physical interpretation of all terms. 
Here Jf (x) is spatial transport of large-scale kinetic en- 
ergy, — PV-u is large-scale pressure dilatation, D({x.) is 
viscous dissipation acting on scales > £, and e^'-'Tx) is 
the energy injected due to external stirring (see [l| for 
details). We have proved in [ij that Di{x.) is negligible 
at scales £ ^ £fj,, where £fj. denotes the dissipation scale. 
We have also shown that mean kinetic energy injection 
can be localized to the largest scales L ^ £ hy proper 
stirring. Over an intermediate scale-range L ^ £ ^ £^^ 
the only relevant terms in eq.Q are inertial processes. 
The subgrid scale (SGS) flux terms are defined as 

nf(x) = -pdjUiT{ui,Uj) (4) 

A,(x) = ^d,Pr{p,u,), (5) 

and act as sinks in ([3]), transferring large-scale kinetic 
energy to scales < £. Deformation work, 11^, is due to 
large-scale strain, Vu^, acting against turbulent stress, 
'p^Te{u,u) ~ p((uu)^ — UfU^), while baropycnal work, 
Af, is due to large-scale pressure-gradient force, ^Pp/pe, 
acting against turbulent mass flux, Ti(p,u) (see [l| for 
a more detailed discus sion of the physics). We employ 
the notation Ti{f, g) = {fg)g — f [g^ for 2"''-order central 
moments of fields /(x), (^(x) 0. 

There are three facts crucial for proving scale locality 
of inter-scale transfer. First is the observation that defor- 
mation work, Ilf, and baropycnal work, A^, represent the 
only two processes capable of direct transfer of kinetic en- 
ergy across scales. Pressure dilatation, — P^V-u^, does 
not contain any modes at scales < £ and does not van- 
ish in the absence of subscale fluctuations. It, therefore, 
cannot participate in transferring kinetic energy directly 
across scales and only contributes to conversion of large- 
scale kinetic energy into internal energy. This observa- 
tion allows us to circumvent analyzing the internal en- 
ergy budget which does not couple to large-scale kinetic 
energy via viscous dynamics, as we have proved in 

The remaining two parts of our proof build upon previ- 
ous studies in incompressible hydrodynamic Q and mag- 
netohydrodynamic [8| turbulence, with some technical 
modifications. The second ingredient we use is the fact 
that SGS kinetic energy flux across £, He + Ae, can be 
written in terms of increments, (5/(x;r) = /(x+r) — /(x), 
for separation distances |r| < £ (or some moderate mul- 
tiple of £) and do not depend on the absolute field /(x). 



Baropycnal work, Af , can be expressed in terms of incre- 
ments by noting that gradient fields and central moments 
are related to increments as 

Vj,^0[Sf {£)/£], f',^0[Sfi£)], 

ri{f,g)^0[Sfi£)Sg{£)] (6) 

where symbol O stands for "same order-of-magnitude 
as," and f'i = f — fe t^i^ fine-scale field. For rigor- 
ous details, see 0, Q- In order to express deformation 
work, n^, in terms of increments, we need the following 
identities which are straightforward to verify: 

Vu = Vu + p-^VT{p, u) - p^'^t{p, u) Vp, (7) 
?(u, u) = r(u, u) + t{p, u, u)/p - t{p, u)r(p, u)/^^. (8) 

We are finally able to express He in terms of increments 
using ([6]) and two additional relations, 

Vr,{f,g) = 0[6f{£) Sg{£)/£], 
Teif,g,h)^0[S.n£) Sg{£) 5h{£)], (9) 

whose rigorous details are in our longer work Q . The re- 
lation of S'^'-order central moments, r (/, g,h) = (fgh)^ — 
fiTi{g, h)-gfTt{f, h)~hfTfif, q)- f ^Qthp, to increments 
is unpublished and due to G. L. Eyink [10|. 

Since Hp and A^ can be expressed in terms of velocity, 
pressure, and density increments, it thus becomes suffi- 
cient to show that these increments themselves are scale- 
local. To establish this, we need the third requirement 
crucial for locality — that scaling properties of structure 
functions of velocity, pressure, and density increments 
are constrained by: 

||du(r)||p ~ Ur,nsAp{r/LYv^ < < 1, (10) 

\\5P{v)\\j, ^ PrmsBp{r/L)<, < 1, (11) 

\\5pir)\\p < prmsCpir/Ly'^, < (12) 

for some dimensionless constants Ap, Bp, and Cp. The 
root-mean-square of a field /(x) is denoted by frms = 
^j2y/2^ where (...) is a space average. Here, the p-th 
power of an Lp-norm |j • \\p = (| ■ |p) is just the tradi- 
tional structure function. We remark that condition 
on the scaling of density increments is only an upper 
bound. It only stipulates that the intensity of density 
fluctuations decays at smaller scales, which is a very mild 
requirement and is readily satisfied in incompressible or 
nearly-incompressible fiows. Heuristically, assumptions 
(fT0ll - (fT2|) characterize the roughness of fields: cr-' < 1 
specifies that the field /(x) is "rough enough", while 
CT-'' > states that /(x) is "smooth enough." 

Under conditions (fT0| - (fT2|) . proving scale-locality of 
the SGS fiux becomes simple and follows directly from 
scale-locality of increments Q. For example, the con- 
tribution to any increment Sf{£) from scales A > ^ is 
represented by 5f^{£). Here, /(x) can denote either ve- 
locity or pressure field. Since the low-pass filtered field 
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(x) is smooth, its increment may be estimated by Tay- 
lor expansion, and ([6]), and (fTO|) or ((TT|) . as 



II^/aWIIp^II^-(v/z 



(13) 



and this is negligible for A 3> ^ as long as ct^ < 1. The 
notation 0{. . . ) denotes a big-O upper bound. On the 
other hand, the contribution to any increment Sf{£) from 
scales S < £ is represented by Sfg{£). Here, /(x) can 
denote either velocity or density field. Since f'g = 0[f{x+ 
6) — f{x)] from scaling conditions ([TUl) and imply 
that 



<2||/^ 



5\\P 



o 



(14) 



and this is negligible for 5 <C ^ as long as > 0. For 
more details and for the careful proofs of these state- 
ments, see 3 and our longer work 

Notice that, unlike for the velocity and pressure fields, 
we do not stipulate that p{x.) be "rough enough." Contri- 
butions to the flux across scale i from the largest density 
scales L ^ £ need not be negligible, yet the flux can still 
be scale-local. The underlying physical reason is simple; 
an energy flux across scale f at a point x will depend on 
the mass in a ball of radius £ around x. Mass is pro- 
portional to average density, ^^(x). in the ball which is 
dominated by large scales: Pf(x) = 0[p^(x)] = 0[prms]- 
Indeed, for incompressible turbulence, the only density 
scale present is a k = Fourier mode, the largest possi- 
ble, and the scale-local SGS flux is directly proportional 
to this density mode. Furthermore, we do not require 
that the pressure field be "smooth enough" even though 
we expect > 0. This is because pressure only ap- 
pears as a large-scale pressure-gradient in with no 
contributions from scales S £. 

The ultimate source of scaling properties pUjl - p^ is 
empirical evidence from experiments, astronomical obser- 
vations, and numerical simulations. For incompressible 
turbulence, which may be viewed as a limiting case of our 
analysis, the scaling of velocity and pressure structure 
functions (fTO|l . (fTT|) has been well-established by a variety 
of independent studies such as those by [ll|, [13 • 
sumption (|12p on density structure functions is trivially 
satisfied for a uniform density field. As for compressible 
turbulence, the available data is also in compelling sup- 
port of our assumptions. Astronomical observations by 
jl3l Il4| show that 2nd-order density structure function 
scales with (Tj == 0.3. Measurements of 2nd-order velocity 
structure functions in molecular clouds 15|, |16[ and solar 
wind [17| also yield < ct^ < 1 for 1 < p < 6. Alongside 
observational evidence, numerical studies of compressible 
turbulence [isl . [Toj report power-law scaling exponents 
well within our required constraints. 



Under an additional assumption concerning the co- 
spectrum of pressure dilatation, which is, albeit rea- 
sonable, not as weak as scaling conditions (fTIH) - (fT^ . 
our proof of a scale-local SGS flux implies a scale- 
local conservative cascade of mean kinetic energy de- 
spite the latter not being an invariant. The require- 
ment on pressure dilatation co-spectrum, E^^{k) = 
Efc-o.5<|k|<fc+o.5-^(k)V-u(-k), is that it decays fast 
enough at large k, 



E G Urms Prms i^L) 



/3>1. (15) 



Here, C is a dimensionless constant and L is an integral 
scale. In the limit of large Reynolds number, assump- 
tion (jlSp implies that mean pressure dilatation, PD{£) = 
(PfV-Uf), converges to a finite constant, 9 = (PV-u), 
and becomes independent of £. In other words, we have 
for wavenumbcr K « £~^, 



lim PD(£) 



lim 



J2 EP^{k)^e. (16) 



We remark that condition (fTS]) is sufficient but not nec- 
essary for the convergence of PD{£) in the limit of £ ^ 0. 
The series X]fc<K E^^{k) can converge with isT — > oo at 



a rate faster than what is implied by assumption (|15p due 
to indcfiniteness in the sign of E^^{k). 

Saturation of mean pressure dilatation in (|16|) reveals 
that its role is to exchange large-scale mean kinetic and 
internal energy over a transitional "conversion" scale- 
range. At smaller scales beyond the conversion range, 
mean kinetic and internal energy budgets statistically de- 
couple. In other words, taking — > first, then £ — > 0, 
steady-state mean kinetic energy budget becomes. 



(H, + A,) = (e™^)-0. 



(17) 



We stress that such a decoupling is statistical and does 
not imply that small scales evolve according to incom- 
pressible dynamics. Small scale compression and rarefac- 
tion can still take place pointwise, however, they yield a 
vanishing contribution to the space-average. 

We denote the largest scale at which such statistical 
decoupling occurs by £c. It may be defined, for instance, 
as 4 = Ek'^"^^^^(k)/Ek-E^^^(k)- Over the ensuing 
scale-range, £c > £^ £p,, net pressure dilatation does not 
play a role, and if, furthermore, (e™-' ) in ([TT]) is localized 
to the largest scales as shown in [ij, then (H^ 4- A^) will 
be a constant, independent of scale £. 

A constant SGS flux implies that mean kinetic energy 
cascades conservatively to smaller scales, despite not be- 
ing an invariant of the governing dynamics. This is one of 
the main conclusions of this Letter. In particular, kinetic 
energy can only reach dissipation scales via the SGS flux, 
H£ H- A^, through a scale- local cascade process. We are 
therefore justified in calling scale-range £c > £ ^ £^ the 
inertial range of compressible turbulence. 
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Needless to say, the scaling of pressure dilatation co- 
spectrum is easily measurable from numerical simula- 
tions. We note that condition (|15p does not require a 
power-law scaling — only that E^^ik) decays at a rate 
faster than ~ It is not at all trivial why one should 
expect PD{£) = {PeV-ue) to converge at small scales. 
How can this be reconciled with the expectation that 
compression, as quantified by V-u, would get more in- 
tense at smaller scales? Indeed, [20| observed numeri- 
cally that (V'u)rms is an increasing function of Reynolds 
number. The key point here is that our assumption (jlSp 
concerns spatially averaged pressure dilatation. It is true 
that V-u(x), being a gradient, derives most of its contri- 
bution from the smallest scales in the flow. Since PV-u is 
not sign-definite, however, major cancellations can occur 
when space-averaging. The situation is very similar to 
helicity co-spectrum in incompressible turbulence, where 
the pointwise vorticity, u;(x) = Vxu, can also become 
unbounded in the limit of infinite Reynolds number. Yet, 
numerical evidence shows that (u-u;) remains finite and 
the helicity co-spectrum decays at a rate ~ with 
n PS 5/3 > 1 (see for example 
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We can offer a physical argument on why PD{1) is 
expected to converge for £ ^ result of cancel- 

lations from space-averaging. The origin of such can- 
cellations can be heuristically explained using decorrela- 
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tion effects very similar to those studied in 
While the pressure in (PV-u) derives most of its con- 
tribution from the largest scales, V-u is dominated by 
the smallest scales in the flow. Therefore, pressure varies 
slowly in space, primarily at scales ~ L, while V-u varies 
much more rapidly, primarily at scales <C i, leading 
to a decorrelation between the two factors. More pre- 
cisely, the pressure Pi in PD{£) may be approximated 
by Pi = 0[Prms] = 0\Pl\ such that 



(PfV-Uf) 



PLV.((u,)^ + (Uf)i)) 

(PiV.u^) + (PL)(V.(u,)i^). 



The first term in the last expression follows from (u^ )^ w 
Ui, while the second term is due to an approximate 
statistical independence between Pl and V-(u£)^ ~ 
Su{£)/£ which varies primarily at much smaller scales 
^ £ L. If there is no transport beyond the domain 
boundaries or if the flow is either statistically homoge- 
neous or isotropic, we get (V-(u£)'^) = 0. The heuristic 
argument finally yields that pressure dilatation, 

PD{£) = (PeV-ue) « (PlV-ul), (18) 

becomes independent of ^, for £ ^ L. Expression pB]) 
corroborates our claim that the primary role of pressure 
dilatation is conversion of large-scale kinetic energy into 
internal energy and does not take part in the cascade 
dynamics beyond a transitional "conversion" scale-range. 

In summary, we conclude that there exists an iner- 
tial range in high Reynolds number compressible turbu- 



lence over which kinetic energy reaches dissipation scales 
through a conservative and scale-local cascade process. 
This precludes the possibility for transfer of kinetic en- 
ergy from the large-scales directly to dissipation scales, 
such as into shocks, at arbitrarily high Reynolds numbers 
as is commonly believed. We make several assumptions 
and predictions which are amenable to empirical scrutiny. 
Our locality results concerning the SGS flux can be veri- 
fled in a manner very similar to what was done in Q and 
[1] . We also invite empirical tests of assumption (fT5|) on 
the scaling of pressure dilatation co-spectrum. Prelimi- 
nary numerical results by (23 | of compressible isotropic 
turbulence indicate that indeed the co-spectrum decays 
at a rate faster than k~^. Verifying ([15]) or (fTH]) under 
a variety of controlled conditions would substantiate the 
idea of statistical decoupling between mean kinetic and 
internal energy budgets. This would have potentially sig- 
nificant implications on devising reduced models of com- 
pressible turbulence. 
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